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ABSTRACT 

We employ three-dimensional state of the art magnetohydrodynamic models of the early solar wind 
and heliosphere and a two-dimensional model for cosmic ray transport to investigate the cosmic ray 
spectrum and flux near the Archean Earth. We assess how sensitive the cosmic ray spectrum is to 
changes in the sunspot placement and magnetic field strength, the large scale dipole magnetic field 
strength, the wind ram pressure, and the Sun's rotation period. Overall, our results confirm earlier 
work that suggested the Archean Earth would have experienced a greatly reduced cosmic ray flux 
than is the case today. The cosmic ray reduction for the early Sun is mainly due to the shorter solar 
rotation period and tighter winding of the Parker spiral, and to the different surface distribution of 
the more active solar magnetic field. These effects lead to a global reduction of the cosmic ray flux at 
1 AU by up to two orders of magnitude or more. Variations in the sunspot magnetic field have more 
effect on the flux than variations in the dipole field component. The wind ram pressure affects the 
cosmic ray flux through its influence on the size of the heliosphere via the pressure balance with the 
ambient interstellar medium. Variations in the interstellar medium pressure experienced by the solar 
system in orbit through Galaxy could lead to order of magnitude changes in the cosmic ray flux at 
Earth on timescales of a few million years. 

Subject headings: Sun: evolution - (Sun:) solar-terrestrial relations - (ISM:) cosmic rays - planets and 
satellites: atmospheres 



1. INTRODUCTION 

During the Archean eon (~ 3.8 — 2.5Gya) when life on 
Earth is thought to have first emerged, the Sun's rotation 
period was 6-15 days, in contrast to its current 27 day 
period (Skumanich 1972). A rotation-activity relation 
is well known from observations of young cool stars, in 
which the faster the rotation the higher the m agnetic ac- 
tivity and related UV- X -ray luminosity (e.g. Pallavicini 



et al. |[i98T Giidel|2007|) . Young active stars are also ob 



served to have magnetic spots concentrated at high lat- 
itudes (i.e. polar regions, Strassmeier||200i I , in contrast 
to the low-latitude spots that characterize the magnetic 
activity of the present day Sun. Since the topology of 
the solar surface magnetic field governs the topology o f 
the interplanetary magnetic field (IMF, Parker 1958b), 
and the topolog y of the solar wind (see summary by 
McComas 2007), the interplanetary environment of the 
young Sun would have be en very different from the one 
at present. Parker (1958a) showed that cosmic ray (CR) 
transport within the solar system, and the flux of ener- 
getic particles that reaches earth, depends on the IMF. 
Consequently, even in the absence of temporal and spa- 
tial variations in Galactic CR production, we can surmise 
that the Archean earth likely experienced a quite differ- 
ent CR exposure to that of the present. 

As an ionizing radiation source, CRs could have played 
a significant role in the origin and evolution of life on 
Earth. Their potential influence includes inducing a 
varied and complex organic chemis try, including pos 



sible production of nucleotides (e.g. Court et al. 2006 
Simakov et al. 2002), causing cellular radiation damage 
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and mutation through direc t and indirect processes (e.g. 
Nelson||2002| |Dartnell|[20lT ) aiding atmo spheric lighten- 
ing initiation (e.g. Gurevich et al. 1999), and (perhaps 
more controversially) affecting clima te through produc- 
tion of cloud condensation nucle i (e.g. Svcnsmark & Frii& 
Christensen|[T997l |Shaviv| [2003l |Wallin;um||20U l| IShaviv 
20051 IKirkby et al. | |2011J )" See aW [Medvedev fc Melott 
(2007n5F an overview. Understanding the way in which 



solar magnetic activity has influence the CR exposure of 
Earth through time is then of interest from a variety of 
both physical and biological perspectives. 

The modulation of galactic CRs in the heliosphere 
is described by the diffusive transport equation of 
Parker ( 1965 ) . The omnidirectional phase space density 



f(v,p, t) at position r, momentum, p, and time, t obeys: 



dt 



V( K V/)-(V + V d )V/ 



VV df 

3 d In p ' 



(1) 



where n stands the anisotropic diffusion tensor with 
different components parallel and perpendicular to the 
Parker-spiral magnetic field (km, and n±, respectively). 
V is the convective solar wind velocity, while de- 
scribes large-scale gradient and curvature drifts due the 
regular gyration of CR particles in the non-uniform field. 
The last term accounts for the energy change due to a 
deceleration in the expanding wind and acceleration at 
t he Terminat i on Shoc k (TS). 



Svensmark ( 2006a 



calculated the CR spectrum at 
1 AU back in time using scaling laws for the solar wind 
and mass loss rate as a function of time, and a scaling law 
for the average solar surface magnetic field as a function 
of rotation rate. Due to the lack of direct observations 
of solar-like winds of cool stars, it is not clear to what 



extent the scaling l aws u sed by jSvensmark (2006a 



are 



valid. Wood et al. ( |2002[ ), for example, caution against 



2 



extrapolating t heir secular wind m ass loss relation that 
was used in the Svensmark ( 2006a ) study to significantly 
higher magnetic ac t ivity levels. 



The Wood et al. (2002) relation is based on a correla- 
tion between X-ray and wind fluxes for a handful of stars, 
and ob s erved secular decline of stellar X-ray luminosity. 
Cohen (2011) argued that stellar mass loss rate should 
not be scaled with stellar X-ray luminosity, since they are 
determined by different components of the stellar mag- 
netic flux (X-rays being associated with closed flux, and 
mass loss with open flux). In short, wind mass loss rates 
for th e early Sun are very uncertain. |Sterenborg et aL| 
(2011) calculated a series of scenarios for the IMF and 
the solar wind structure for the young Sun. Using a mag- 
netohydrodynamic wind model and a surface magnetic 
field consistent with observations of young active stars, 
they found a wind mass loss rate not more than 10 times 
that of the current Su n. This result is consistent w ith 
the theoretical study of Holzwarth k. Jardine| fl2007|) . 
Here , we extend the work done by |ijterenborg et al. 



(2011 ), and build on their wind and IMF model solutions 



to investigate the CR spectrum near the Archean Earth 
using a two-dimensional model for transport of CR in 
the heliosphere. In particular, we focus on the effect 
of latitudinal variations of solar active regions and the 
enhanced solar rotation. We do not assume anything 
about the effect of the TS and the heliopause, or the 
location of the solar system within the Galaxy. In other 
words, we assume the same incoming CR distribution for 
each possible solar magnetic field topology. 

In Section [2] we describe our experimental approach 
and the assumptions being made in both modeling the 
solar wind and the IMF of the young Sun, and the cal- 
culation of CR transport. We describe the results in 
Section [3] and discuss their meaning in Section [4] 

2. MODEL DESCRIPTION 

Our calculation of the CR flux near the Archean Earth 
involves building a model of the solar wind and the IMF 
for the young Sun, which was done by [Stcrcnborg 'et al.| 
(2011 ), and the calculation of the CR transport from the 
edge of the solar system to the Earth based on the output 
from the solar wind model. The two steps are as follows. 

2.1. Solar Wind and the Interplanetary Magnetic Field 

The transport of CR depends mostly on the topology 
of the magnetic field in the heliosphere (i.e., the IMF). 
The IMF is governed by the surface distribution of the 
solar magnetic field, which is carried away by the solar 
wind, and w ound up by the solar rotation. |Sterenborg| 
et al.j (|2011[) have used t he BATS-R-US Magnetohydro- 



1999 Toth et al. 



dynamic s~[MHD) code (Powell et al. 
2005| Cohen et al.||2007 l to calculate the solar wind and 
the IMF conditions tor the young Sun, based on a set 



young 

of scenarios for the young Sun's surface field topology, 
as well as for different rotation rates. For each solution, 
a synoptic map of the solar surface field (magnetogram) 
was used to calculate a steady-state, three-dimensional 
solution for the solar corona and the solar wind. 

For the young Sun, a magnetogram of the current Sun 
during solar maximum was used as the basis for a ref- 
erence case. The magnetogram was manipulated and 
split into two components — a weak dipole component, 



denoted by a Dipole Factor (DF), and a strong spot 
components, denoted by a Spot Factor (SF). In order 
to resemble observations of young active stars, the active 
regions in the map were shifted by 30 and 60 degrees 
towards the poles, so that they appeared at mid-high 
and high latitudes. The case of 60 degrees shift should 
represent the solar magnetic field at the begining of the 
Archean , while the 30 degrees shift should represent 
the field to wards the end of that period (e.g., Toner & 



Gray 



Figure [T] shows the magnetograms and the 
three-dimensional magnetic field structure for the refer- 
ence case of the current Sun and for the reference young 
Sun cases with the spots manipulated and shifted pole- 
ward. Other test cases involved scaling the magnitude of 
either magnetogram components (DF or SF) by a factor 
of 10. 

Here, we use these solutions to calculate the transport 
of CR from the edge of the solar system to 1AU. Since 
we are interested in the average state of the young Sun's 
IMF, longitudinal variations are not important. This of- 
fers some computational expedience. We found the wind 
solutions are fully developed and super- Alfvenic at a ra- 
dial distance of 20Rq, such that a much more computa- 
tionally expensive fully-detailed calculation out to 1AU 
is not necessary for understanding the longitudinally- 
averaged wind at at Earth. We therefore extract the 
physical parameters from the model at 20i?©, averaged 
over longitude, and extrapolate the model to the TS. We 
have computed models for rotation periods of 2, 4.6, 10, 
and 26 days to probe the effects of rapid rotation and 
winding of the Parker spiral. 

The beginning of the Archean, about 3.8 Gya, corre- 
sponds to solar ages comparable to open clusters such 
as M37 (550 Myr) and the Hyades (625 Myr), when the 
sola r rotation period wo uld have been about 4-6 days 
(e.g. Meibom et al.|20lT ), while throughout the Archean, 
the solar rotation pe riod was increasing to 10-15 days 
( |Telleschi et "al][2005l ) . 

2.2. Cosmic Ray Transport 

The modulation of CR is simulated in an updated ver- 
sion of th e quasi-steady-state , two dimensional numerical 
model by Jokipii et al. ( 1993). An average uniform radial 
solar wind with the value of the average wind speed is 
taken from eac h of the different young Sun models de- 
scribed in Sect. [O The TS is placed at 90 AU for the 
current Sun, and at distances scaling with the square 
root of the average ram-pressure for other cases (95 AU, 
110 AU, and 120 AU for polar spots, enhanced DF, and 
enhanced SF, respectively). We perform additional cal- 
culation of all cases with the TS remain fixed at 90 AU 
for reference. 

We neglect azimuthal gradients, and use the az- 
imuthally averaged values of the magnetic field for the 
diffusion coefficients, Ky and k±, as well for the drift 
velocity, V^. Though this is not a strictly valid as- 
sumption, it serves as a good approximation for our 
purposes to estimate the omnidirectional CR flux. For 
the diffusion coefficients, we adopt the simplest scaling 
law, assuming that both k» and k± are proportional to 
the particle speed times the gyroradius, i.e. they are 
scaled inversely with the large-scale field, B. We take 
K|| = 18.1P(P/B) AU 2 /day : and k± = 0.05Ky. Here, 
/8 = v/c is the particle speed as a fraction of the speed 
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of light, P = pcjZe is the particle rigidity, which is es- 
sentially the momentum per unit charge, and B is the 
magnetic field strength in units of nT. We also add a 
transverse field component impeding easy penetration 
through the weak spiral field throu gh the polar regions 
according to Jokipii & Kota (1989). 

The omnidirectional density, /, relates to the particle 
flux, Jt, through Jt — Anfp 2 , and Jt is set at t he outer 
boundary to the in terstellar spectrum given by [Webber 
& Lockwoodl pOOll as: 



21. IT- 2 - 8 



1. + 5.85T" 1 - 22 + 1.18T" 2 - 54 



(2) 



with T being the particle kinetic energy. 

One omission of the present work is the treatment of 
CR transport through the heliopause and the TS itself, 
since we are interested in the CR transport within the 
solar system. 

3. RESULTS 

The results for the CR energy spectrum at 1 AU for 
the simulated test cases are shown in Figure [2] Each 
plot shows the simulated CR spectrum for the different 
solar rotation periods. The plots also show the effect of 
variations in the location of the TS for the cases rep- 
resenting the young Sun. The result of the spectrum 
for the current Sun with the current rotation per iod are 
consis t ent wi th those presented in Figure 1 in |Svens- 
mark (2006a I ( presenting both model calculation and 
data taken from |McDonald et al. (2001)), with the flux 
peak around 1 — 2 [particles m~ z s~ L sr~ 1 MeV~ 1 ] lo- 
cated around 300 - 500 MeV. 

The most dramatic change, which is clearly seen in all 
panels, is the change in the spectrum due to the enhanced 
solar rotation. The faster rotation winds up the field car- 
ried by the wind and increases the azimuthal component 
of the IMF. The enhanced and compressed Parker spiral 
prevents CR penetration to 1 AU much more efficiently 
than the Parker spiral of the current Sun. Even for the 
relatively low activity of the present day Sun, the en- 
hanced rotation reduces the flux at the peak of the spec- 
trum by more than an order of magnitude. It also shifts 
the spectrum peak by a factor of 2 towards higher ener- 
gies from about 500 MeV to 1 GeV. Figure [3] shows the 
equatorial magnetic field for the cases with solar rota- 
tion periods of 4.6 and 26 days, along with a conceptual 
display of spirals with different rotation rates. While the 
difference in the field tangling is notable in the simulation 
box extended up to 24i? , it becomes larger at greater 
heliospheric distances. 

The other effect that is clearly seen from the results is 
that of the sunspot distribution, and of the magnitude 
of the different field components. By only moving the 
spots by 60 degrees towards the pole, a slight — 30-50% 
or so — but global decrease in the CR flux is obtained. 
This change is due to the increase in the polar field as 
the spots occupy high-latitude regions. Since it is easier 
for CR to penetrate through polar regions, where the 
field is more radial and relatively weaker, an increase in 
the polar field leads to a reduction of the CR flux. The 
shift by 30 degrees does not seem to change the CR flux 
at all, probably since the shifted spots do not interact 
and modify the polar field as in the 60 degree shift case. 



Figure [4] shows the results for the cases with enhanced 
dipole and spot components of the magnetic field. When 
the dipole component of the reference young Sun polar 
spot model is enhanced by a factor of 10 (a polar field 
of about 50 G), the spectra for the slow to fast rotation 
periods are suppressed by factors of about 2 and 4-5, re- 
spectively, with only very modest changes in the spectral 
peaks. For the case where the polar spot component is 
enhanced by a factor of 10 (field strengths on small spa- 
tial scales of few kG), both slow and fast rotation spectra 
are decreased much more dramatically — by almost an or- 
der of magnitude. In fact, all CR with energies less than 
100 MeV are completely eliminated from the spectrum 
for the cases with 4.6 and 2 days rotation period. 

The location of the TS does not have a large influence 
on the spectrum, except for the case of enhanced polar 
dipole/spots where it is moved from 90 AU (the current 
Sun) to 120 AU. The TS being further out causes the 
CR to travel longer distances from the edge of the solar 
system to 1 AU. This results in a larger CR loss on the 
way and an overall reduction in the CR flux that reaches 
the inner solar system. The change in the CR spectrum 
for the other models is smaller partly because their wind 
solutions are more similar in terms of wind momentum 
to that of the present day active Sun case and the TS 
therefore not so far out. For the polar spots model, CR 
coming from high latitudes are also attenuated more be- 
cause of the significantly enhanced polar field. 

4. DISCUSSION AND CONCLUSIONS 

Our simulations demonstrate how the different inter- 
planetary environment during the early Archean would 
have resulted in a quite different CR spectrum and 
greatly reduced flux near the Earth than we experience 
today. What we have done here is to isolate the different 
aspects of the magnetic activity of the Sun to determine 
to what extent each affects the resulting CR spectrum. 
The CR reduction we find is mainly due to the enhanced 
solar rotation and the different surface distribution of 
the more active solar magnetic field. These effects lead 
to a global reduction of CR flux at 1 AU by up to two 
orders of magnitude or more, and to elimination of the 
low-energy (< lOOMeV) flux in the most extreme 2 day 
rotation period scenario. The results are more sensitive 
to the assumed strength of polar spots than to the large- 
scale dipole field. 

4.0.1. Comparison with previous work 



Shaviv (2003) estimated that a stronger solar wind 



from a more active young Sun would effect a reduction 
in the CR flux reaching Earth that could have influenced 
the climate. The results presented here illustrate that a 
stronger wind, giving a larger heliosphere and TS shock 
distance is not the main effect; instead the more rapid so- 
lar rotation that winds the Parker spiral is the dominant 
screening mechanism. 

It is useful t o compare our resul ts with those presented 
in Figure 1 of Svensmark ( 2006a[ ); we have adopted the 
same units for displaying the cosmic ray energy spec- 
trum as that study. We generally confirm their findings, 
though with some differences. In their curve for 3.8 bil- 
lion years ago, almost the entire CR energy spectrum is 
eliminated except for energies above about 10 GeV where 
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the flux is reduced one hundredfold. This is similar to 
our last case in Figure R1 with a strong polar (spot) field 
and a rotation period of 2 days — more extreme, at least 
in terms of rotation, than we believe the Sun at that age 
would have been. As noted in Sect. [2] at the begining of 
the Archean (the Eoarchean era), the solar rotation pe- 
riod would have been 4 -6 days, which seem s somewhat 
slow er than adopted by Svensmark] ( |2006a ) . 



In Svensmark (2006a 



the effect of the magnetic field 
distribution was treated by simply scaling the average 
surface field with the rotation rate, with a resulting field 
strength as a function of time, t, of B(t) = B t~ 6 , with t 
in units of 4.6 Gyr. At I Gyr, this formula predicts a field 
2.5 times that of today — a fairly modest difference. Our 
simulations enable us to separate the effects of rotation 
rate, field distribution, field strength, and the location of 
the TS. Therefore, they provide intermediate solutions 
as well. In some of our scenarios, the CR flux does not 
become negligible, especially in t he lower energy ran ge. 



This contrasts with the results of Svensmark ( 2006a ) 
4.1. Surface Field of the Archean Sun 

The sensitivity of our results to the assumed polar field 
strength — comparing the panels in Figure[2]and[4] — raises 
the question of what polar spot field strength is appro- 
priate for the young Sun? 

The polar regions of stars are difficult to probe in de- 
tail with Doppler imaging techniques. These methods 
employ rotationally-modulated wavelength shifts in Zee- 
man signatures to m ap the surface field structure ( |Donati| 
fc Landstreetj200"9 ) . Polar regions exhibit only small pro- 



jected line-ot-sight velocity changes through a rotation 
cycle (or none at all at the pole itself), and consequently 
the resolution of the mapping technique is much more 
limited. In particular, a polar region with mixed polarity 
field, as characterizes our young Sun reference model, will 
unlikely register as a strongly magnetic region because 
the opposite polarity Zeeman signatures will cancel out. 
Moreover, the dark polar spots commonly seen on active 
stars contribute little to the full-disk spectrum and can 



Marsden et al. 


2006a 


Marsden et al. 


2006b 



Indeed, parsden et al. (2006a), 



and Jetters & Donati J 200 



found 
hundred G 

circling a dark polar spots within which no significant 
field was detected on the 30-50 Myr old, solar analogs 
HD 171488 (P ~ 1.3d) and HR 1817 (P ~ l.Od). At 
somewhat later spectral types that might not be as rele- 
vant to the solar case, the well-studied very active zero- 
age main-sequ ence K0 d warf AB Doradus (age ~ 75 Myr; 



Janson et al. 



e.g. 



20071 shows some mixed polarity re- 
gion s at 50 deg latitudes with field strengths up to 500 G 

(e.g. |Hussain et aL|2007[ |Cohen et aL|2010|) . 

Larg er scale fields are more easily resolved. |Petit et ah] 
(2008) found a mean large scale field of 51 ± 6G for the 
solar analog HD 190771, which has a rotation period of 
8.8d. This is an order of magnitude larger than that of 
the Sun, with one third the rotation period. 

It seems likely, then, that both the large scale dipo- 
lar and smaller scale spot fields on the young Eoarchean 
Sun were significantly stronger than seen on the Sun of 
today. The model with a 4 day rotation period and spot 
and fields enhanced by a factor of 10 is then probably not 
too far from the truth. While this model produces peak 



field amplitudes of more than a kG, the stellar obser- 
vations are limited in spatial resolution and would tend 
to smear and average out field with fine structure. As 
the Sun aged and was rotationally braked by the solar 
wind through the Archean, the spot and large scale field 
strengths would have declined and shifted towards lower 
latitudes. 

4.1.1. Broader Relevance 
The relevance of CR for modulation of the Earth's 



climate — first raised as a possibility by Ney ( 1959 ) — has 
not yet been placed on a firm footing! Early evidence 
suggest i ng a C R-climate connection was summarised by 
|Shaviv| ( |2003[) . The general mechanism is thought to be 
the seeding of cloud condensation nuclei by CR ionization 
events in the lower atmosphere that results in a higher 
albedo and cooler average temperatures. While several 



studies have refuted the con nection (e.g. Erlykin et al. 
|2009[ |Pierce fc Ada ms 2009 , and references therein), re 



cent experimental support for the bas ic mechanism has 
been presented by Kirkby et al. (2011). 

Over all, our results co nfirm the findi ngs of |Svensmark| 
( 2006a ) and conjecture of Shaviv ( 2003 ) that the Archean 
Earth would have experienced a greatly reduced CR flux 
than is the case today. If CR do play a significant role 
in cloud seeding, it seems plausible that the paleoclimate 
of early Earth could have been affected by reduced cloud 
formation. 

4.1.2. The Termination Shock Distance 

While solar rotation plays the dominant role in sec- 
ular CR modulation, the timescalc for its evolution is 



hundreds of Myr. |Svensmark (2006b I further suggested 
that the CR variations would occur as the Sun passed 
through regions of enhanced star formation during its or- 
bit through the Galaxy. One additional CR modulation 
effect probed by our models is the variation in interstellar 
medium (ISM) pressure in which the solar system finds 
itself. The heliopause contracts and expands according 
to the balance between the ISM pressure and the solar 
wind ram pressure. The radius of the heliopause is then 
roughly proportional to the square root of the local ISM 
pressure. 

The ISM is largely in pressure equilibrium, with the 
pressure balance shared between gas thermal and dy- 
namic pressures, magnetic fields and cosmi c rays (see, 
e.g., the reviews of |Ferriere||2001| |Cox||200"5| . However, 
the gas density and temperature distribution is highly in- 
homogeneous, and, although these components are also 
in approximate pressure balance, significant variations in 
pressure are thought to exist. Self-gravitating molecular 
clouds can have pres s ures s everal times that of the am- 



bient medium. Cox (2005) assesss the spiral arm and 



inter-arm region average pressure differences to be a fac- 
tor of 2 and p erhaps more. In a study of the local bubble 
and environs, |Lallement et al. (2003) notes that the gas 
pressure in "local fluff" in which the Sun is currently 
located is between 5 and 10 times less than the pres- 
sure deduced from soft X-ray background measurements 
and that magnetic pressure making up the differences 
is "hardly co mpatible with the observat ions" . However, 
later work by Koutroumpa et al. (20091 has shown that 
much of the soft X-ray background is from foreground 
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heliospheric emission, therefore, the pressure of the lo- 
cal bubble surrounding the solar system is still highly 
uncertain. 

The current picture of the ISM then indicates the Sun 
will pass through pressure differences of factors of per- 
haps up to ten. For a relative velocity of the order of 
25 km s" 1 — approximat ely the current solar velo city rel- 
ative to the local ISM (Lallement et al. 2005) — and a 
scale size of significant pressure inhomogeneities of, say, 
50 pc, the timescale for CF flux variations would be of 
the order of a few Myr. 

A factor of two in pressure corresponds to a change in 
TS radius of 40%, and our models indicate a CR flux 
change of a factor of 2 or so. A factor of ten pressure 
change would lead to CR flux changes of more than an 
order of magnitude. CR fluxes would be expected to 
be higher during spiral arm passages, where the aver- 
age pressure is higher, but could also spike when passing 
through molecular clouds whose higher pressures could 
greatly compress the heliosphere. The mechanism of 
heliospheric compression would augment (and perhaps 
dominate at times) the s piral arm pas sage increase in 
the ambient CR flux that |Shaviv ( 2005| ) an d |Svensmark| 
(2006a I have suggested could have le ad to ~ 250 Myr cy- 



on the heliospheric size and structure using an MHD 
model, which includes the effect of neutral atoms. Based 
on their model cases for the heliosphere, they calculated 
the change in CR intensity at 1 AU and found that the 
flux can differ by an order of magnitude, depending on 
the density and ionization state of the ISM region in 
which the solar system finds itself. Compression of the 
heliosphere due to higher ISM pressure generally leads 
to larger CR fluxes at Earth in their models, consistent 
with the d irection of the effec t modeled much more sim- 
ply here. (|Fields et al. 2008) found that the passage of 



the blast from a supernova explosion lOpc away would 
shrink the heliopause down to about 1 AU. Such "as- 
trospheric collapse" has been estimated to occur with a 
frequency of 1-10 Gyr -1 d ue to passage throug h very 
dense interstellar clouds by (S mith fc Scalo|2009 1 . 



cies in the climate and biodiversity. Miiller et al. ( 2006 1 
have studied the effect of the change in ISM pressure 
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Fig. 1. — The observed surface distribution of the photospheric radial magnetic field during solar maximum is shown on the top left, 
while the manipulated magnetogram with the spots shifted by 30 and 60 degrees toward the poles is shown on the top middle and right 
panels, respectvely. The three-dimensional magnetic field topologies corresponding to these three surface field distributions (displayed on 
each sphere) are shown on the bottom. 
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Fig. 2. — CR energy spectrum for solar rotations of 26d (current rotation, blue), lOd (green), 4.6d (magenta), and 2d (red). The local 
ISM spectrum is shown in black. Plots are for the current Sun (left), spots shifted towards the pole by 30 degrees (middle), and spots 
shifted towards the pole by 60 degrees (right). Solid lines represent the spectrum with the TS scaled with the solar wind dynamic pressure, 
while dashed lines represent the spectrum for the TS fixed at 90 AU. 



Fig. 3. — Left: conceptual display of different spirals with slow (blue) and fast (red) rotation rates as a function of distance. Right: the 
simulated equatorial magnetic field near the Sun (up to 24Rq) for the case with solar rotation period of 4.6 days (red) and 26 days (blue) 
displayed from the top. The difference in tangling increases with distance from the Sun. 




